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The presence of macroscopic phase separation into superconducting and spin-glass-like 
magnetic phases in LaFeAs(Oi_. E F. I ;) is demonstrated by muon spin rotation measurement 
in a sample near the phase boundary (x — 0.06). Both magnetism and superconductivity 
develop simultaneously below a common critical temperature, T m ~ T c ~ 18 K. This re- 
markable accordance strongly suggests that the electronic correlations leading to these two 
competing ground states share a common origin. 
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The recent discovery of the oxypnictide superconductor LaFeAsOi_ x F x (LFAO-F) with 
a critical temperature (T c ) of 26 K 1 and the successful revelation of much increased T c upon 
the substitution of La for other rare-earth elements (such as Sm, leading to ~ 43 K 2 ) and 
the application of pressure for LFAO-F (~ 43 K 3 ) have triggered broad interest in the mech- 
anism yielding a relatively high T c in this new class of compounds. They have a layered 
structure like high-T c cuprates, where the dopant and conducting layers are so separated that 
the doped carriers (electrons introduced by the substitution of O 2- with F~ in the La2C>2 
layers) move within the layers consisting of strongly bonded Fe and As atoms. They exhibit 
another qualitative similarity to cuprates in that superconductivity occurs upon carrier doping 
of pristine compounds that exhibit magnetism. 4 Some preliminary results of the muon spin 
rotation/relaxation (/^SR) experiment on a variety of oxypnictide superconductors showed 
that the superfluid density n s falls on the empirical line on the n s vs T c diagram observed for 
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the underdoped cuprates, 5,6 from which possibility of the common mechanism of supercon- 
ductivity is argued between oxypnictides and cuprates. 

However, in terms of the doping phase diagram, there are certain differences between 
these two systems, e.g., (i) T c (> for 0.4 < x < 0.12) does not vary much with x 1 as 
in cuprates known as "bell-shaped" and (ii) the magnetic (spin density wave, SDW) phase 
shares a boundary with the superconducting phase near x ~ 0.04. 1,7 The insensitivity of T c 
to x is reasonably understood from the conventional BCS theory where condensation energy 
is predicted to be independent of carrier concentration. The close relationship of magnetism 
and superconductivity suggests that a detailed investigation of how these two phases coexist 
(and compete) near the phase boundary will provide important clues to elucidating the paring 
mechanism. Among various techniques, ^SR has a great advantage in that it can be applied 
in systems consisting of spatially inhomogeneous multiple phases, providing information on 
respective phases according to their fractional yield. Our /xSR measurement in the LFAO-F 
sample with x = 0.06 (T c ~ 18 K) reveals that these two phases indeed coexist in the form of 
macroscopic phase separation, and more interestingly, that a spin glass-like magnetic phase 
develops in conjunction with superconductivity in the paramagnetic phase. This accordance 
strongly suggests a common origin of the electronic correlation leading to these two competing 
phases. 

Although the oxypnictide with rare-earth (R) substitution i?FeAsOi_a;F x exhibits higher 
T c than that of LFAO-F, strong random fields from rare-earth ions preclude a detailed study 
of the ground state using sensitive magnetic probes like /uSR. Therefore, we chose the original 
LFAO-F system for our ^xSR study. The target concentration of LaFeAsOi^F^ is set near the 
phase boundary, x = 0.06, for which a polycrystalline sample was synthesized by solid state 
reaction. The detailed procedure for sample preparation is described in an earlier report. 1 
The sample was confirmed to be mostly of single phase using X-ray diffraction analysis. Of 
two possible impurity phases, namely, LaOF and FeAs, only the latter exhibits a magnetic 
(helical) order with Tjy — 77 K. 8 As shown in Fig. 1, magnetic susceptibility exhibits no trace 
of FeAs phase or local magnetic impurities except below ~ 50 K where a small upturn is 
observed. The susceptibility at a lower field [shown in Fig. 3(a)] provides evidence of bulk 
superconductivity with T c ~ 18 K from the onset of diamagnetism. Conventional /xSR mea- 
surement was performed using the LAMPF spectrometer installed on the M15 beamline of 
TRIUMF, Canada. During the measurement under a zero field (ZF), residual magnetic field 
at the sample position was reduced below 10~ 6 T with the initial muon spin direction parallel 
to the muon beam direction [P M (0) || z]. For longitudinal field (LF) measurement, a magnetic 
field was applied parallel to P^(0). Time-dependent muon polarization [G z (t) = z ■ P^,(t)] was 
monitored by measuring decay-positron asymmetry along the i-axis. Transverse field (TF) 
condition was realized by rotating the initial muon polarization so that P M (0) || x, where the 
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Fig. 1. (Color online) Magnetic susceptibility of LaFeAsOi-^F^ with x = 0.06 for the sample used 
for /xSR measurement. Inset shows a reduced view of the region below 35 K. 

asymmetry was monitored along the x-axis to obtain G x (t) = x ■ P^(t). All the measurements 
under a magnetic field were made by cooling the sample to the target temperature after the 
field equilibrated. 

ZF-//SR is the most sensitive technique for examining magnetism in any form, where the 
development of local magnetic moments leads to either the spontaneous oscillation (for long- 
range order) or exponential damping (inhomogeneous local magnetism) of G z {t). Figure 2 
shows examples of ZF-/xSR time spectra collected at 2 and 30 K. The spectrum at 30 K 
(> T c ) exhibits a Gaussian- like relaxation due to weak random local fields from nuclear 
magnetic moments, indicating that the entire sample is in the paramagnetic state. Meanwhile, 
the spectrum at 2 K is split into two components, one that exhibits a steep relaxation and 
the other that remains to show Gaussian-like relaxation. This indicates that there is a finite 
fraction of implanted muons that sense hyperfine fields from local electronic moments. The 
absence of oscillatory signal implies that the hyperfine field is highly inhomogeneous, so that 
the local magnetism is characterized by strong randomness or spin fluctuation. The fractional 
yield of the component showing steep relaxation is as large as 25% (see below), which is 
hardly attributed to impurity and therefore implies that the sample exhibits a macroscopic 
phase separation into two phases. 

The magnitude of the hyperfine field and that of spin fluctuation are evaluated by observ- 
ing the response of the ^SR spectrum to a longitudinal field (LF). It is shown in Fig. 2 that 
the relaxation in the paramagnetic component is quenched by applying a weak magnetic field 
(LF=5 mT), which is perfectly explained by the suppression of static nuclear dipolar fields 
(< 10° mT). Meanwhile, the faster relaxation (seen for < t < 1 /xs) due to the magnetic 
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Fig. 2. (Color online) /iSR time spectra observed in LaFeAs(Oi_ a; F :E ) with x — 0.06 at 2 K under 
a longitudinal field (LF), a zero field (ZF), and a transverse field (TF), and that under ZF at 30 
K. The spectrum under TF is plotted on a rotating reference frame of 6.78 MHz to extract the 
envelop function. 

phase is recovered only gradually over a field range of lO 1 ^ 2 mT, and there still remains a 
slow relaxation even at the highest field of 60 mT. This residual depolarization under LF is 
a clear sign that local spins are slowly fluctuating, leading to the spin-lattice relaxation of 
Pu(t). Such quasi-two-step relaxation is also observed in dilute spin-glass systems, 9 which is 
understood as a distribution of spin correlation time. A detailed analysis is made considering 
that these two components coming from the magnetic phase (see below). 

Under a transverse field, implanted muons experience an inhomogeneity of the field [5 z (r)] 
due to flux line lattice formation below T c that leads to relaxation, in addition to those ob- 
served under a zero field. The TF-/iSR time spectrum in Fig. 2 (envelop part of the oscillation) 
obtained under a field of 50 mT exhibits complete depolarization at 2 K, indicating that the 
entire volume of the paramagnetic phase falls into the superconducting state. The rapidly re- 
laxing component observed under ZF is also visible (although the coarse binning of the spectra 
for extracting the envelop part makes it slightly obscure), indicating that the corresponding 
part of the sample remains magnetic. 

Considering the presence of the magnetic phase besides the paramagnetic ^supercon- 
ducting below T c ) phase, we take special precaution to analyze both TF and ZF/LF /j,SR 
spectra in a consistent manner. For the determination of physical parameters describing the 
behavior of signals from the magnetic phase, we first analyze ZF/LF spectra at 2 K using the 
X-square minimization fit with the relaxation function 

3 

G z (t) = H + Y j w i exp(-A i t)] -G KT (5 N : t), (1) 

8=2 
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where Gkt(^n : t) is the Kubo-Toyabe relaxation function for describing the Gaussian damp- 
ing due to random local fields from nuclear moments (with <5n being the linewidth), 10 w\ is 
the fractional yield for the paramagnetic phase, w 2 and W3 are those for the magnetic phase 
E w i = -0 w ith A2 and A3 being the corresponding relaxation rate described by the Redfield 
relation 

A, = -fe (< = 2,3), (2) 

i ' n 

where w M = 7^-Hlf, 7^ is the muon gyromagnetic ratio (= 2ir x 135.53 MHz/T), Hjj? is the 
longitudinal field, 5 2 and #3 are the means of the hyperfine fields exerting on muons from 
local electronic moments, and v 2 and z/3 are the fluctuation rates of the hyperfine field. The 
solid curves in Fig. 2 show the result of analysis where all the spectra at different fields 
(only ZF and LF=5mT are shown here) are fitted simultaneously using eqs. (1) and (2) with 
common parameter values (except cj m that is fixed to the respective value for each spectrum), 
which show excellent agreement with all the spectra. The deduced parameters are as follows: 
wi = 0.754(9), w 2 = 0.165(9), w 3 = 0.081(4), S 2 = 0.71(5) /is -1 , 5 3 = 3.9(3) /is -1 , v 2 = 1.7(2) 
//s -1 , and v 3 = 4(1) //s -1 . Although the depolarization in the magnetic phase is approximately 
represented by two components with different hyperfine couplings (Si), the fluctuation rates 
(ui) are close to each other (10 7 s _1 at 2 K), suggesting that the randomness is primarily due 
to the distribution in the size of local moments (or in their distances to muons). Since no 
impurity phase with a fraction as large as 25% is detected by X-ray diffraction analysis, it is 
concluded that this magnetic phase is intrinsic. 

In the analysis of temperature-dependent TF spectra, we used the relaxation function 

G x (t) = exp(-i^t 2 )[t i ;iexp(-<5 s 2 t 2 )cos(27r/ s t + (/>) 

+ (w 2 + w 3 )exp(-A m t)cos(27r/ m t + 0)], (3) 

where Wi and <5n are fixed to the values obtained by analyzing ZF/LF-/iSR spectra. The first 
component in the above equation represents the contribution of flux line lattice formation in 
the superconducting phase, where S s corresponds to the linewidth a s = \/2S s = 7 M ((£?(r) — 
Bq) 2 ) 1 / 2 [with B being the mean B(r)], while the second term represents the relaxation in 
the magnetic phase. Here, the relaxation rate for the latter is represented by a single value 
A m (instead of A2,3), as it turns out that the two components observed under LF are hardly 
discernible in TF-^SR spectra. [This does not affect the result of the analysis, because the 
amplitude is fixed to w 2 + W3 so that A m may represent a mean ~ (w 2 A 2 + W3A3) / (w 2 + W3).] 
The fit analysis using the above form indicates that all the spectra are perfectly reproduced 
while the partial asymmetry is fixed to the value determined from ZF-//SR spectra. This 
strengthens the presumption that the paramagnetic phase becomes superconducting below 
T c . The result of analysis is summarized in Fig. 3, together with the result of dc magnetization 
measured in the sample from the same batch as that used for ^SR. 
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It is interesting to note in Fig. 3(b) that, although the central frequency in the supercon- 
ducting phase (/ s ) does not show much change below T c ~ 18 K probably owing to a large 
magnetic penetration depth (it is indeed large, see below), that in the magnetic phase (/ m ) 
exhibits a clear shift in the negative direction below T m ~ T c . The magnitude of the shift is 
as large as ~ 1% and thus is readily identified despite a relatively low external field of 50 mT. 
As shown in Fig. 3(c), the relaxation rate in the magnetic phase (A m ) also develops below 
T c in accordance with the frequency shift, demonstrating that a spin-glass- like magnetism 
sets in below T c . Here, we note that the development of magnetic phase is already evident 
in the ZF/LF-/iSR spectra, and results are fully consistent with each other. The onset of 
superconductivity below T c is also confirmed by an increase in <5 S , as observed in Fig. 3(c). 
This remarkable accordance of onset temperature between magnetism and superconductiv- 
ity strongly suggests that there is an intrinsic relationship between the superconducting and 
magnetic phases that leads to a common characteristic temperature. 

The temperature dependence of a s in Fig. 3(d) is compared with theoretical predictions 
for a variety of models with different order parameters. The weak-coupling BCS model (s- 
wave, single gap) apparently fails to reproduce the present data, as they exhibit a tendency 
to vary with temperature over the region T/T c < 0.4. Although a fit using a two-gap model 11 
shown by a dotted line seems to exhibit reasonable agreement with data, the deduced gap 
parameters (2Aj//cb7c) are largely inconsistent with the prediction of the weak-coupling BCS 
model (see Table I). These observations suggest that the superconducting order parameter is 
not described by a s-wave symmetry with a single gap. When a power law, a s = [1 — (T/T c )^] , 
is used in fitting the data, we obtain a curve shown by the broken line in Fig. 3(d) with an 
exponent (3 ~ 2. This is in good agreement with the case of <i-wave symmetry at the dirty 
limit. 

In the limit of extreme type II superconductors [i.e., A/£ S> 1, where A is the effective 
London penetration depth and £ = \J <&q/ (2itH C 2) is the Ginzburg-Landau coherence length, 
$o is the flux quantum, and H C 2 is the upper critical field], a s is determined by A using 
the relation 12 ct s / 7m = 2.74 x 10" 2 (1 - h) [l + 3.9(1 - h) 2 ] 1/2 $ A~ 2 , where h = H TF /H c2 



Tabic I. Parameters for defining the lines in Fig. 3(d). 



Two-gap 


Power law 


T C (K) 


18.0(5) 


UK) 


18.9(4) 


(7(0) GUS- 1 ) 


0.291(5) 


(7(0) G^s- 1 ) 


0.291(4) 


w 


0.73(6) 





1.7(1) 


2Ai/fe B T c 


2.6(3) 






2A 2 /fe B T c 


1.1(3) 
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Fig. 3. (Color online) Temperature dependence of dc magnetic susceptibility measured at lmT (a), 
and that of physical parameters deduced by analyzing TF-/^SR spectra in superconducting (/ s , 
5 S ) and magnetic (/ m , A m ) phases (b-c), and of <t s (= V2S S ) proportional to superfluid density 
(d). Curves in (d) are fits by models described in the text. 



and Htf is the magnitude of external field. From a s extrapolated to T = and taking 
H c2 ~ 50 T (ref. 13 ), we obtain A=595(3) nm. Because of the large anisotropy expected from the 
layered structure of this compound, A in the polycrystalline sample would be predominantly 
determined by in-plane penetration depth (A a b). Using the formula A = 1.31A a b for such 
a situation, 14 we obtain A a b=454(2) nm. This value coincides with that expected from the 
aforementioned empirical relation between A~ b 2 superconductors. 6 ' 15 However, this may not 
be uniquely attributed to the superfluid density because A depends not only n s but also on 
the effective mass, a s oc A~ 2 = n s e 2 /m*c 2 . 
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Finally, we discuss the feature of the spin glass-like phase. Assuming that the local mo- 
ments are those of off-stoichiometric iron atoms with a moment size close to that in the SDW 
phase (~ 0.36/Ub 4 )) the mean distance between muon and iron moments in the relevant phase 
is estimated to be ~ 0.5 nm from an average of 8{. Given the unit cell size (a = 0.403 nm, 
c = 0.874 nm 1 ), this would mean that more than a quarter of iron atoms in the magnetic 
phase (i.e., ~ 7% of the entire sample) should serve as a source of local moments. It is unlikely 
that such a significant fraction of iron atoms remains as impurities in the present sample. 

It might also be noteworthy that there is an anomaly near T m 2 — 12 K in the susceptibility 
[the onset of ZFC/FC hysteresis in Fig. 1 and a steplike kink in Fig. 3(a)]. This seems to be 
in accordance with the onset of a steeper increase in A m below T m 2, suggesting a change in 
magnetic correlation. 

/[/SR studies of LFAO-F have been made by a number of groups. According to those 
preliminary studies, no clear sign of magnetism is observed in the sample over relevant doping 
concentrations, except for a weak relaxation observed far below T c for x = 0.05 and 0.075 
in ZF-fiSK spectra and an unidentified additional relaxation observed in TF-/iSR spectra for 
x = 0.075. 5-7 This led us to recall the sensitivity to chemical stoichiometry in the emergence 
of the spin glass-like A-phase observed near the boundary between the antiferromagnetic 
and superconducting phases in CeCu2Si2- 16 In addition to the vl-phase, the present LFAO-F 
system exhibits a closer similarity to this classical heavy-fermion superconductor such as the 
phase diagram against pressure/doping. 17 Further study of the dependence of fractional yield 
for the magnetic phase with varying x (in small steps near the phase boundary) would provide 
further insight into the true nature of these phases and the mechanism of superconductivity 
itself that is working behind the coexistence/competition. 

In summary, it has been revealed by our //SR experiment that superconducting and mag- 
netic phases coexist in LaFeAs(Oo.94Fo.o6) with x = 0.06. These two phases simultaneously 
develop just below T c , strongly suggesting an intimate and intrinsic relationship between 
these two phases. The result of TF-/uSR measurement suggests that the superconductivity 
of LaFeAs(Oo.94Fo.o6) cannot be explained by the conventional weak-BCS model (single gap, 
s-wave) . 
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experiment. This work was partially supported by the KEK-MSL Inter-University Program 
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